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The parametric amplifier with and without the pumping fluctuations of coupling function 
is considered when the fields are initially prepared in coherent light. The pumping fluctu- 
ations are assumed to be normally distributed with time-dependent variance. The effects 
of antibunching and anticorrelation of photons on the photon distribution, correlation be- 
tween modes and factorial moments are demonstrated. A possible enhancement of photon 
antibunching for certain values of initial mean photon numbers is shown and discussed. We 
have shown also that new states (called modified squeezed vacuum states or even thermal 
states) can be generated from such an interaction. Further, we have demonstrated that the 
sum photon-number distribution can exhibit collapses and revivals in the photon-number 
domain somewhat similar to those known in the Jaynes-Cummings model. 

PACS numbers: 42.50Dv,42.60.Gd 



I. INTRODUCTION 



Quadratic parametric processes have been considered since 1961, when a quantum- mechanical 
model of the parametric amplification and frequency conversion, treating the pump field classically 
and neglecting the losses, had been proposed in [1]. Among these processes parametric amplifier 
takes a considerable interest since degenerate and nondegenerate parametric amplifiers can perfectly 
provide single-mode and two-mode squeezing (squeezed light has less noise than a coherent light 
in one of the field quadratures and can exhibit a number of features having no classical analogue). 
Further, parametric amplifier has been employed in experiments. For example, the fourth-order 
interference effects arise when pairs of photons produced in parametric amplifier are injected into 
Michelson interferometers Q|. The second-order interference is observed in the superposition of 
signal photons from two coherently pumped parametric amplifiers when the paths of the idler 
photons are aligned j^]. 
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It is well-known that the interaction of light with matter changes the statistical properties of 
light. The nature of the change depends on the kind of interaction process and on the statistical 
properties of the incident electromagnetic field. An investigation of the statistical properties of the 
nondegenerate parametric amplifier (without losses) has been performed using different techniques, 
such as Heisenberg technique [4] and Schrodinger technique ja|. The main attention in these articles 
has been paid to the single mode case, i.e. only one mode (signal or idler) is detected. The 



6], 



where the variance of the photon number 



anticorrelation in this model is an interesting effect 
is less than the average photon number and the photocounting distribution becomes narrower 
than the corresponding Poisson distribution for a coherent state with the same mean photon 
number. These nonclassical phenomena are established for certain values of the phase of incident 
fields. On the other hand, the quantum statistical properties of the parametric amplification 
and generation including losses in the medium have been studied by means of the generalized 
Fokker-Planck equation for the distribution related to the antinormal ordering [7] showing that the 
normal counting generating function decomposes into the product of two functions which are the 
generating functions for the superposition of coherent and chaotic fields and that the outgoing field 
exhibits the anticorrelation effect. Further, it is convenient mentioning that the photon counting 
distribution together with the waiting-time distribution techniques have been used to analyze the 
nonclassical effects for light beams generated by the lossy degenerate optical parametric amplifier 
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12. HQ modelled 



Finally, we refer to parametric processes with stochastic pumping 
by means of a stochastic differential equations. One can find the total photon number, as a function 
of time, of frequency up-conversion process has a turning point and the antibunching effect is 
predicted regardless of the mean photon numbers available in the process. Also the effect of 
pumping fluctuations are considered in 



12, 



13, 



14|], where the observed sub-Poissonian behaviour 



is progressively degraded as the level of fluctuations in the pumping radiation increases. This 
means that the pump fluctuations play a similar role as reservoir, where the interaction with the 
environment leads to the destruction of the nonclassical effect, i.e. to the coherent energy extraction 
from the quantum mechanical system. 

In the present paper we study the quantum properties of the parametric amplifier with the 
classical pumping by conside ring the pump fluctuations in the system. More illustratively, the 



lossless effective Hamiltonian 
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14] controlling the system is 
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H = hwia\ai + hu>2a-2a2 — hg[aid,2 exp(iu>t — i<j>) + h.c], (1) 

where a,j,j = 1,2, are the annihilation operators designated to the signal and idler mode, respec- 
tively; u)j are the natural frequencies of oscillations of the uncoupled modes with uj = u\ + 0J2 
(assuming the resonance frequency condition holds); <f) is the initial phase of the pump and h.c. 
is the Hermitian conjugate; g is the real coupling constant (phase matching is assumed) and we 
assume g = go + egi, that is the coupling constant consists of two parts: go is the deterministic 
part and g\ is the correction part (i.e. |<?o| >> ksi|)> an d e is a real fluctuating quantity obeying 
time-dependent Gaussian distribution described by 



P(e,t) 



y/2ira(t) 



cxp 



2a(t) 



(2) 



where the standard deviation is a(t) = <ro[l — exp(— fj,t)], ao is the maximum standard deviation 
and fi is a finite real parameter. This is the difference from what has been done earlier 12. [13I. 



where the standard deviation is time independent. In this way pump fluctuations are taken into 
account. Indeed, time-dependent standard deviation can lead to the enhancement of nonclassical 
effects compared to the time-independent case, as we shall see later. For any observable F(e,t), 
the average processs will be performed as follows: 



((F(e,t))) = / (F(e,t)) { P(s,t)dE, (3) 



where the two angle brackets means that there are two average operations, one is the average 
(expectation value) over the initial states of the field (..)f and the other is the average over the 
stochastic parameter e using the distribution ([2]) in the sense described on the right-hand side of 
([3]). Such an averaging process will be frequently used in this article. 

We should stress here that the model will be treated generally for two cases with and without 
the pump fluctuations where the comparison is instructive. Also we give new insight into the 
possibility to generate nonclassical light based on the change of intensity of initial incident fields 
which has not been considered earlier [3, 3, [3, Q] • 

For completeness, the solution of the Heisenberg's equations for the Hamiltonian (JTJ) is well 
known 
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&i(t) = exp(— iu>it)[ai(0) cosh(gt) + ia 2 (0) sinh(gt) exp(icp)], 

(4) 

<h(t) = exp(— iu>2t)[a,2(0) cosh(gt) + ia[(0) smh(gt) exp(i^)]. 

The small term \e\g\ « \go\ is considered as a small perturbation and the standard deviation 
a(t) is assumed as a slowly varying quantity compared to the variations in (|4|). 

This article is organized as follows: In section 2 the basic equations and relations such as 
anticorrelation functions, reduced factorial moments and sum photon- number distribution are given 
when the modes are initially prepared in coherent light. In section 3 discussion of the results is 
performed and the conclusions are summed in section 4. 

II. BASIC RELATIONS AND EQUATIONS 

As it is well known light fields for which the P representation is a well-behaved distribution 
cannot exhibit nonclassical features [if]) . P-function for the single mode as output from para- 
metric amplifier ([I]) is well defined, because in the parametric systems the evolution broadens the 
single mode P distribution due to the spontaneous pump photon decay [13]. This also increases 



181 ] and reflects that there is no 



the radius of the Wigner contour compared to the initial one 
(single mode) nonclassical behaviour. For this reason we focus our attention on the behaviour of 
the compound modes when the modes are initially prepared in coherent states. 

Now, cross-correlation may be used to describe anticorrelation between modes, which may cause 
that the variance of the photon number is less than the average of the photon number giving then 
antibunching; it can be measured by detecting single modes separately by two photodetectors and 
correlating their outputs. Cross-correlation between mode 1 and mode 2 is controlled by 



(: AW X AW 2 :) = ((: al(t)ai(t)aT(t)a 2 (t) :)) - ((a{(t)a 1 (t))(aT(t)a 2 (t))), (5) 

where Wj = a^(t)a,j(t),j = 1,2, and : : denotes the normally ordered operator, i.e. creation 
operators aj are to the left of annihilation operators aj; the double angle brackets in the right- 
hand side has the same meaning as before. 



an nn 



For input coherent light (| 01,02)) equation © together with (JH) reduce to [71. 1 12. 1 la. 

(: AW1AW2 :) = ^(2|ai| 2 + 2|a 2 | 2 + l)sinh 2 (2 5 t) 
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«i I \o.2 1 sinh(4gi) sin ip 



-(2|ai| 2 + 2|a 2 | 2 + 1) {cosh(4 5o t) exp[8a(t) ff 2 t 2 ] - l} 



|a;i I |ck2 I sinh(4<7o^) exp[8cr(i)<? 2 i 2 ] simp 



(6) 



(7) 



where atj = \aj\exp(i(f>j),j = 1,2 ( \atj\ and cftj are the amplitude and the phase of the initial 
coherent jth mode (signal or idler)), and ip = <ft\ + 4>2 — (ft represents the difference between the 
phase-sum of the initial signal and idler modes and the initial phase of the pump, which is the 
quantity determining the phase mismatch. In © and in the following we have performed the 
average with respect to the initial coherent states and the average over stochastic pumping field, 
which is indicated by the angle brackets (in ([7]) the integration over e has been carried out using 
the distribution ([2]) after prescription ([3])). 

We proceed by writing down the formulae of the reduced factorial moments and sum photon- 
number distribution for compound modes; the details of the derivation for these relations can be 
found in [lj|, respectively as 



(8) 



k\fk-l\fl 



k k—l I k—l\—l 2 

= £££ £ 

Z=0 /i=0 Z 2 =0 h=0 

k-h- l 2 \ (-l) k+l3 ~ l - h (k - l)W.h l ^h 2 2 
l 3 J hlhl 

X (l)k-h-l2 exp [ 2go t(l -h- h) + 2a(t)g 2 1 t 2 (l - h - h) 2 } 



P(n,t) 



exp 



Aj A 2 1 

JT+XT) (1+A 2 )J 



(1 + Ai)(l + A 2 



<n - 



A 2 q + Ai y-« 



1=0 
xL n _/[- 



(n -/)!/! 
Ai 



(1 + A 2 
A2 , 



Ai(l + Ai) J lL A 2 (l + A 2 
2^i=o 2^h=o 2^i 2 =o 2^i 3 =o 2^/4=0 I h ) \i 2 ) \ 



v^' v^»-/j x (n— 1\ ( I \ (n—h—h ^ x +l 2 +l 



(9) 

(10) 
(11) 



(-l)"+'3+'4-'-'i( n+ ; 1+ ; 2+ ; 4+ i)!fe' 1 i/ i ^ 
hH 2 \hl{n+h+l 2 +iy. 
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where 



(12) 



Ai = \ [exp(-2 5 i) - 1] , A 2 = \ [exp(2gt) - 1] , A l = h 1 exp{-2gt) , 
A 2 = h 2 exp(2gt), hi = |(|ai| 2 + \a 2 \ 2 + 2|ai||a 2 | sba.ip), 
h 2 = ^(|ai| 2 + \a 2 \ 2 — 2|ai||a 2 | sim/>), 

p /.\ _ A(h 1 -h 2 )tgi cxp(-2gpt) r r.\ _ r,, _ 2Qi-fe 2 ) 

[l+cxp(-2 So t)] 2 ' L 2\ L )- zn l l+cxp(-2g t) ' 

r(i, ij) = Ti(t) + 2^it(Zi + Z 2 + h + k + 1), 



r (t, Zj) = r 2 (t) + 2 5o t(/i + z 2 + h + z 4 + 1) 

and L n (.) is the Laguerre polynomial of order n; as before the angle brackets in (|8|) and (jlOp 
represent the average over the stochastic pumping variable e. For simplicity we have dropped the 
time t from the coefficients Aj and \j. The quantities Aj and Xj play the role of the mean numbers 
of coherent photons and mean numbers of chaotic photons, respectively The transition from (|8j) 
to ([9]) as well as from (|10[) to (jlip has been done by carrying out the integration over the stochastic 
variable. It is reasonable mentioning that the expressions Q and are different from those in 
Qj [l3 ] , where the average over the initial complex field amplitudes is taken there and consequently 
additional summations have appeared. When performing this integration the following expressions 
have been used in the treatment of the exponential factor for the sum photon-number distribution: 



+ M-=2hi- 2(/ii - h 2 )[l + exp(-2 5 t)] 



-i 



l+Ai 1 1+A 2 



2/ii - 2(h! - h 2 ) EZoi- 1 ) 1 exp(-2Z 50 t)[l - 2 9l tle + 0(e 2 )} (13) 



2hi - 2(h x - h 2 



l+exp(-2 90 i) + [l+exp(-2 90 t)] 2 

where in these relations g is replaced by its value, and 0(e 2 ) includes terms of the power equal 
or higher than e 2 , which have been neglected in the last step since e is small. In spite of this 
the neglection is physically sensible and necessary for performing the average over the parameter 
e, it leads to a problem that the distribution (jlip does not exist as an ordinary well behaved 
function. Consequently additional treatments may be used to avoid this difficulty. Therefore, if 



2g\tle cxp(— 2got) 
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the propagation time through the crystal, which is also the interaction time, is sufficiently short 
together with the fact that the contribution of the correction part (the term involving g\) must 
be small compared with the main contribution, some modification of the correction term can be 
performed. For example we can approximate g 2 t 2 l1 by g 2 t 2 {l\ + I2 + Z3 + 1)^4 in the exponential 
term. The choice of the factor (l\ + I2 + Z3 + 1) is to compensate the large values of Z4 as much 
as possible since lj,j = 1,2,3 are the indices of the summations (cf. (jllf) ) . It is important to 
point out for a short time interaction (as well as finite values of g\ and <7o) and finite values of 
I4 that the contribution of r|(i) is considerably small compared with that of To(t). The validity 
of this approximation can be checked numerically. Now the sum photon-number distribution (jlip 
becomes well behaved function, as it should be, and having the form 

p(„ +\ _ V" V 1 " 1 V* s^n-h-h (n-l\ ( l \ (n-h-l 2 \ (-lr+'s-'-'ife'/fe^i+^+i 
r\'hl) — 1^1=0 l^h=0 Z^J 2 =0 2^/3=0 Wi )\l 2 )\ l a ) hU 2 l 

-(n+h+h+2) 



x < 1 + exp 



2a(t)g^(l + r -ff + 3(h+l 2 + k + l))-2tg \y' n 1 2 (14) 



x exp {-2t 50 (Z 1 + l 2 + h + 1) - r 2 (t) + ^(t)[ri(i) + 2t gi (h + l 2 + h + I)] 2 } • 

One can easily check that when gi — > 0, all equations for fluctuating pumping reduce to the 
nonfluctuating ones. 

In the following we investigate the behaviour of the system on the basis of the results of the 
present section. 



III. DISCUSSION OF THE RESULTS 



The discussion of the results of the model under consideration will be performed through two 
parts. In the first part we give an investigation for both the anticorrelation between modes and 
the reduced factorial moments. However, the second part is devoted to the analysis of the sum 
photon-number distribution. Also in this part we develop new states and discuss some of its 
properties. 



A. Anticorrelations and reduced factorial moments 



We start our discussion in this part by analyzing the anticorrelations between modes in the 
system, which can show some interesting features. For usual parametric amplifier (i.e. g = go), 
one can easily check that the extreme values of © are established at if) = ±7r/2, i.e. the maximum 
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FIG. 1: The normalized correlation of fluctuations between modes 

[K(t) = (: AW 1 AW 2 :)/{Wi){W 2 ) - 1], against \a\ for ip = ir/2,t = 0.5, g = I, fi = = 0.6, for 

the time-dependent fluctuating pumping with ctq=0.1 (solid curve), 0.5 (long-dashed curve) and for the 
time-independent fluctuating pumping with <7o=0.1 (short-dashed curve), 0.5 (circle-centered curve). The 
star-centered curve is corresponding to that of the usual parametric amplifier and the straight line is the 
bound of anticorrelation. 

value is at ip = —tt/2 and minimum value at ip = ir/2. Careful examination shows that it is always 
positive when ip = —tt/2, i.e. the behaviour is always classical. The maximal nonclassical effect 
can appear when ip = tt/2. In other words, the anticorrelation is observable as a result of coupling 
of modes. Further, considering that a\ = ot\ = a, t > 0, and ip = tt/2, one can easily obtain that 
anticorrelation occurs if the relation © is negative and this provides 

\a\ > ^exp(g t)y / smh(2g t); (15) 

if the equality holds this expression means that coherent light can be recovered. This relation 
gives the range of \a\ over which the anticorrelation can be obtained. The right hand side of (|15f) 
is monotonically increasing function in the course of time and then the nonclassical effect can 
be realized for a certain time interval after switching on the interaction. It can be shown that 
the maximum value of anticorrelation between modes is — | exp(4go£) sinh 2 (2(7oi)j which occurs at 
\a\ = \ exp(2<7oi) smh(4got). Also the maximum value of correlation establishes at \a\ = 0. Indeed, 
all these informations are clear in Figs. 1,2, where the normalized forms of both correlation 
between modes and reduced factorial moments have been plotted against \a\, respectively. From 
these figures one can see that the anticorrelation reflects itself in the moments of reduced factorial 
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FIG. 2: The reduced factorial moments ((W k )/(W) k - 1) against |a| when k = 5 for ip = n/2, t = 0.5, g = 
= = 0.4, for the time-dependent fluctuating pumping with <7o=0.1 (solid curve), 0.3 (long-dashed 
curve) and for the time-independent fluctuating pumping with er =0.1 (short-dashed curve), 0.3 (circle- 
centered curve). The star-centered curve is corresponding to that of the usual parametric amplifier and the 
straight line is the bound of antibunching. 

moment. In general, the behaviour of these quantities is smooth where, e.g. for reduced factorial 
moment, all curves start from the bunching values going to the antibunching ones as the initial 
mean photon number increases until they reach a maximum, then they start again to decrease the 
antibunching values and eventually become stable for a large domain of \a\. The most effective 
antibunching can arise for the usual parametric amplifier as expected (see star-centered curves). 
When the pumping field fluctuations are considered, the degradation of the nonclassical effect 
occurs (see Figs. 1 and 2). From these figures, the comparison between the solid and short-dashed 
curves as well as between the long-dashed and circle-centered curves shows that the nonclassical 
effects for the standard deviation in time-dependent case are more pronounced than those for the 
time-independent case. Further, one can see also that as the value of increases, the amount of 
antibunching decreases. 

Now we would like to conclude this part by giving an example when the interaction time is 
short which is relevant for practical situation of propagating waves in nonlinear media. We discuss 
this point only for the normalized correlation between modes where similar arguments can be 
given for the reduced factorial moments. This has been displayed in Fig. 3 for the shown values 
of the parameters. We note when the standard deviation is time-dependent, the behaviour of 
the normalized correlation of fluctuations between modes is the same for all chosen values and 
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FIG. 3: The normalized correlation of fluctuations between modes for ip — tt/2, \a.j\ 2 = 2, /j, = 1, go = 10 7 s 1 , 
for the time-dependent fluctuating pumping (solid curve) and for the time-independent fluctuating pumping. 



7„-i 



1, bell-centered 



Short-dashed curve is for g\ — 10 s~ ,ao = 10, long-dashed curve for g\ = 10 s 
curve for g\ — 10 7 s^ 1 ,<jq = 10. The solid curve is for the usual parametric amplifier as well as for the 
time-dependent fluctuating pumping. 



coincides with the behaviour of the usual parametric amplifier, which gives the most effective 
anticorrelation values (see solid curve in Fig. 3). This seems to be related to the short interaction 
time. Also we noted that the anticorrelation nature becomes less pronounced as go (or interaction 
time t) increases, where the amplification of radiation can occur and this is connected with an 



increase in chaotic noise [20 



Indeed, the regime of anticorrelation or antibunching is related 



to the damping of the field t]j Qj- ^ n conclusion one can control the intensity of initial fields, 
phases, and interaction time in such a way that we can get maximum antibunching effect. This 
phenomenon can be observed by detecting both the modes 1 and 2 simultaneously beyond the 
nonlinear medium. It is reasonable to mention that the sum of mode noises of W = W\ + W2 
can be less than modulus of negative values of 2(: AW1AW2 :), i.e. antibunching can occur if 
(: {AW) 2 :) < ((: AVF1AW2 :) < 0), which we have also verified. Such situation cannot occur 
in classical physics and this leads to the nonexistence of the Glauber-Sudarshan quasidistribution. 
Such quiet light may be applied in optical signal processing and optical communications. 
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B. Sum photon-number distribution 

Firstly at ay =0, j = 1,2, for usual parametric amplifier, equation (|10p for the sum photon- 
number distribution reduces to 

P(,,t) = '- ta ?^f (16) 
cosh (gt) ^ 

This expression shows that only even terms survive and this is nonclassical effect. Moreover, such 
photon-number distribution indicates that the system collapses to produce the state 

!*(*)> = ^sk^EZo[^H9t)) l eMilx)\2l), 

(17) 

= d^J EM,[tanh(^)exp(i X )]i[l + (-l)']|l>, 

where x is a relevant phase. In fact, state (|17p is quite different from both squeezed vacuum states 
21] and single-mode of two-mode squeezed vacuum states [jjj. These two states are determined 



respectively by 



E^-ifSpo, (18) 

y/cosh[gt) 2 /! 



1 oo 

PL® = cosh 2 (gt) E[ ta nh(gt)] 2 ^)(/|, (19) 

where p\(t) is the reduced density operator representing the single-mode case of two- mode squeezed 
vacuum states. 

The comparison of the three states (|17p -(|19 p is instructive. One can find that states (|17p are 
completely different from states (|19p . where the latter is the state usually associated with chaotic 
Bose-Einstein statistics |jjj], i.e. it cannot exhibit nonclassical effects. However, the former can 
exhibit nonclassical effects due to the superposition principle. Actually, states (|17p can be written 
in the form of even thermal states by tracing over the phase through the relation 

(20) 

= iir££o(^) 2 Wz|, 

where n = sinh 2 (gi); in this case n represents the mean number of thermal photons and p et h means 
the density matrix of the even thermal states. On the other hand, states (|17p and (|18|) are different 
in their structures but they can have similar nonclassical effects. For instance, both states can have 
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FIG. 4: The photon-number distribution of squeezed vacuum states (dashed-curve) and of modified squeezed 
vacuum states (solid curve) against the number of photons n for t — 0.5, . go = 2, gi = 0. 



pairwise oscillation in the photon-number distribution but with different amplitudes, also they can 
display quadratures squeezing rather than sub-Poissonian statistics (which can be easily checked). 
So that we shall call state (|17p as a modified squeezed vacuum state (or even thermal state). In 
Fig. 4 we have plotted the photon-number distribution for squeezed vacuum and modified squeezed 
vacuum states for comparison. From this figure it is clear that the oscillations in the photon-number 
distribution of state (|17p are more pronounced than those of (I18p . These oscillations are purely 
quantum effect without classical analogue. The origin of these pairwise oscillations in (|18p is 
two-photon nature of the s quee ze operator and they have been interpreted by means of the 



interference in phase space 



22 ] or in the framework of the generalized superposition of coherent 



fields and negative quantum noise 



231 ] . Nevertheless, for the state (|17p the origin of the oscillatory 



behaviour in the photon-number distribution is the quantum interference between the components 
of the state (see second line of (fl~7|) ) in phase space. In fact, such a type of interference is quite 



different from the interference of the Schrodinger cat states [24 ]. To make this point clear we 
may analyze the behaviour of Wigner (W) function for modified squeezed vacuum states since 
this function is informative and sensitive to the interference in phase space. Applying standard 
technique, the PF-function of (|17p is 



W(z) = —1^ {exp [-gt - 2|z| 2 e-^] + exp [gt - 2\z\^] } , (21) 



where z = x + iy. 
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is, 




FIG. 5: The W-function of modified squeezed vacuum against Rez = i, Imz = y for the same situation as 
Fig. 4. 



This relation shows that the TV-function cannot exhibit negative values and its contour in 
phase space is symmetric in contrast with that of squeezed states showing noise ellipse. Further, it 
contains two terms, the first one is arising from the statistical mixture and the second represents the 
interference between the components of the state (117p . Consequently the behaviour of TV-function 
here depends on the competition between these two terms where both have a Gaussian shape with 
center at the origin, however, the first term is always broader than the second one. This is clear in 
Fig. 5 where TV-function given by (|2ip is plotted. In fact, such interference is responsible for the 
nonclassical effects of the modified squeezed vacuum states. Further, similar behaviour has also 
been seen for even binomial states 



251 ]. which contains a finite number of Fock basis states. Finally, 
it should be stressed here that states (fT7|) cannot be produced by a unitary operator as the squeezed 
vacuum state. This is so becuase there is no Hermitian Hamiltonian from which the state (|17p can 
be generated via a canonical transformation. To be more specific, substituting |2Z) = ^t|2Z) into 
the right-hand side of (|17p . one can find that the summation cannot be performed to obtain the 
required unitary operator and this is in contrast with (|19p . Further, the second line of (|17p shows 
that this state represents a superposition of two like-phase-coherent states and consequently 
it cannot be obtained from any unitary operation [^j. 

Now let us return to the sum photon-number distribution when a>j ^ 0, j = 1,2. Actually 
we noted that it is in a rough agreement with the normalized correlation of fluctuations between 
modes. For instance, for the usual parametric amplifier and when tfi = — vr/2, t > 0, P(n, t) almost 
exhibits single-peak structure broader than that of the corresponding coherent light with the same 
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n n 

FIG. 6: The sum photon-number distribution against the number of photon n for a) usual parametric 
amplifier when ip = n/2,t = 0.5, 51 = 0, fi = 1, (go, |oy|)= (1,2.5884) (short-dashed curve), (1.5,0.87) (long- 
dashed curve), the solid curve is for the corresponding Poissonian distribution with the same (W(t)); b) 
time-dependent fluctuating pumping when V = tt/2, ( = 3x 10~ 4 , g 1 = 0.06, a Q = 0.02, (g , \a 3 ■ |) = (10 4 , 2). 

mean photon numbers. However, when tp = ir/2, P(n,t) can exhibit single-peak or oscillatory 
behaviour according to the relation between the incident mean photon numbers, coupling constant 
and the interaction time (see Figs. 6a and b for shown values of the parameters). In Fig. 6a \a\ = 
I ex.p(2got) sinh(Agnt) , he. when the maximum antibunching effects occur in normalized correlation 
function between modes. One can observe that the distribution for this case (short-dashed curve) 
is narrower than that of the corresponding Poissonian distribution showing nonclassical effects. 
Similar behaviours have been seen for the values of the parameters satisfying inequality (|15p . On 
the other hand, when [a| = \ exp(got)y / smh(2got) — 0.25, i.e. the model exhibits always bunching 
as it is clear from (]15p . In this case the photon- number distribution can exhibit nonclassical 
oscillations (long-dashed curve), however, it is broader than that for Poissonian light (which is 
not included in the figure). The sum photon- number distributions of time-dependent and time- 
independent fluctuating pump have a similar behaviour for short interaction time (see Fig. 6b). 
From this figure there are evident the large scale macroscopic oscillations in the behaviour of 
P(n,t). To be more specific, P(n,t) exhibits collapses and revivals similar to those familiar in 
the Jaynes-Cummings model (JCM) [28]; the former is in the photon number domain rather than 
the time domain. In fact, observations on the collapses and revivals of Rabi oscillations were 
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reported 



29)]. Further, collapses and revivals in the JCM are showing the granular structure of the 



photon-number distribution, however, here they indicate the establishment of strong correlation 
between the signal and idler modes in the nonlinear crystal. Indeed, similar phenomenon has been 

□ n 

seen for the photon-number distribution of single-mode [30] and two-mode [31[ squeezed coherent 



states with complex squeeze and displacement parameters. For the two-mode case the phenomenon 
has been realized for the diagonal, and joint photon-number distribution having a fixed number of 
photons in one of the modes, however, the case displayed here (Fig. 6b) for the sum photon-number 



distribution is much richer than for the above distributions 



3l|. 



IV. CONCLUSIONS 



As was shown earlier [7|, |l2|, LL3[] one can get antibunching from the system described by ([I]) 
at a short time and for specific choice of the parameters. Also such an effect can be obtained 

Q 

from parametric amplification process with quantum pumping [6J regardless of the initial phase 
condition. In this case the source of this effect is the quantization of the pumping. In this 
article we have analyzed the role of the initial mean photon numbers on generated nonclassical 
light (antibunched light) which has not been considered earlier [?], Q, 13, 14] so that the overall 



quantum efficiency for the effect can be enlarged. In other words, the degree of antibunching can be 
controlled depending on the initial mean photon numbers. Further, we have shown that considering 
time-dependent pump fluctuations, they decrease the rate of degradation of the nonclassical effects 
compared with the time-independent case. This is evident from the behaviour of the reduced 
factorial moments and the normalized correlation of fluctuations between modes. The behaviour 
of the sum photon-number distribution is in agreement with the behaviour of reduced factorial 
moments and normalized correlation of fluctuations between modes showing antibunching, however, 
in some situations they can display nonclassical oscillations including their collapses and revivals. 

We have shown also that we can generate a modified squeezed vacuum states (or even thermal 
states) when the two-mode vacuum states evolve in the present interaction, thus detecting the sum 
photon-number distribution for outgoing light using photodetectors. These states have somewhat 
similar behaviours as squeezed vacuum states, however, the origin of the nonclassical effects for 
these states is in the quantum interference between components of the state. The photon-number 
distribution for such states exhibits more pronounced oscillations than those of squeezed vacuum 
states for the same chosen parameters. 



16 



Acknowledgement 

J. P. and F. A. A. E-0. aknowledge the partial support from the Projects VS96028, LN00A015 
and Research Project CEZ: J 14/98 of Czech Ministry of Education and from the Project 
202/00/0142 of Czech Grant Agency. One of us (M. S. A.) is grateful for the financial support 
from the Project Math 1418/19 of the Research Centre, College of Science, King Saud University. 

References 



[1] W. H. Louiscll and A. Yariv, Phys. Rev. 124 (1961) 1646. 

[2] P. G. Kwiat, W. A. Vareka, C. K. Hong, H. Nathel and R. Y. Chiao, Phys. Rev. A 41 (1990) 2910; Z. 

Y. Ou, X. Y. Zou, L. J. Wang and L. Mandcl, Phys. Rev. Lett. 65 (1990) 321. 

[3] X. Y. Zou, L. J. Wang and L. Mandel, Phys. Rev. Lett. 67 (1991) 318. 

[4] B. R. Mollow and R. J. Glauber, Phys. Rev. 160 (1967) 1076; 1097. 

[5] L. Mista, Czech J. Phys. B 19 (1969) 443. 

[6] L. Mista and J. Pefina, Acta Phys. Pol. A 52(3) (1977) 425. 

[7] L. Mista and J. Pefina, Czech J. Phys. B 28 (1978) 392. 

[8] R. Vyas and S. Singh, Phys. Rev. A 40 (1989) 5147; A. B. Dodson and R. Vyas, Phys. Rev. A 47 (1993) 
3396. 

[9] W. S. Kryszewski and J. Chrostowski, J. Phys. A 10 (1977) 261. 

[10] S. K. Srinivasan and S. Udayabaskaran, Opt. Act. 26 (1979) 1535. 

[11] W. J. Miclniczuk, Opt. Act. 26 (1979) 1115. 

[12] V. Pcfinova, Opt. Act. 28 (1981) 747. 

[13] V. Pcfinova and J. Pefina, Opt. Act. 28 (1981) 769. 

[14] J. Pefina: Quantum Statistics of Linear and Nonlinear Optical Phenomena, 2nd edition (Kluwer, 

Dordrecht 1991). 

[15] R. J. Glauber, Phys. Rev. 130 (1963) 2529. 

[16] E. Schrodinger, Naturwiss. 14 (1926) 664. 

[17] A. Bandilla and H.-H. Ritze, Opt. Commun. 34 (1980) 190. 

[18] S. M. Barnctt and P. L. Knight, J. Opt. Soc. Am. B 2 (1985) 467. 

[19] S. M. Barnett and P. L. Knight, J. Mod. Opt. 34 (1987) 841. 

[20] L. Mista, V. Pcfinova, J. Pefina and Z. Braunerova, Act. Phys. Pol. A 51 (1977) 739. 

[21] H. P. Yuen, Phys. Rev. A 13 (1976) 2226. 



17 



[22] W. Schleich and J. A. Wheeler, Nature 326 (1987) 574; W. Schleich and J. A. Wheeler, J. Opt. Soc. 

Am. B 4 (1987) 1715; W. Schleich, D. F. Walls and J. A. Wheeler, Phys. Rev. A 38 (1988) 1177. 
[23] J. Pefina and J. Bajer, Phys. Rev. A 41 (1990) 516; J. Pefina, Z. Hradil and B. Jurco, Quantum Optics 

and Fundamentals of Physics (Kluwcr, Dordrecht 1994). 
[24] E. Schrodinger, Nature 23 (1935) 44. 

[25] M. S. Abdalla, M. H. Mahran and A.-S. F. Obada, J. Mod. Opt. 41 (1994) 1889; A.-S. F. Obada, M. 

H. Mahran, F. A. A. El-Orany and M. S. Abdalla, Int. J. Th. Phys. 35 (1996) 1393. 
[26] J. H. Shapiro, S. R. Shepard and N. W. Wong, Phys. Rev. Lett. 62 (1991) 2377; J. H. Shapiro and S. 

R. Shepard, Phys. Rev. A 43 (1991) 3795. 
[27] K. Wodkiewicz, P. L. Knight, S. J. Buckle and S. M. Barnctt, Phys. Rev. A 35 (1987) 2567. 
[28] E. T. Jaynes and F. Cummings, Proc. IEEE 51 (1963) 89; N. H. Ebcrly, N. B. Norozhny and J. J. 

Sanchez-Mondragon, Phys. Rev. Lett. 44 (1980) 1323. 
[29] G. Rempe, H. Walther and N. Klein, Phys. Rev. Lett. 58 (1987) 353 
[30] B. Dutta, N. Mukunda, R. Simon and A. Subramaniam, J. Opt. Soc. B 10 (1993) 253. 
[31] M. Selvadoray, M. S. Kumar and R. Simon, Phys. Rev. A 49 (1994) 4957. 



